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ABSTRACT 



We analyze PLANET photometric observations of the caustic-crossing binary-lens microlensing event, 
EROS BLG-2000-5, and find that modeling the observed light curve requires incorporation of the microlens parallax 
and the binary orbital motion. The projected Einstein radius (fg = 3.61 ± 0.11 AU) is derived from the measurement 
of the microlens parallax, and we are also able to infer the angular Einstein radius (#e = 1.38 ± 0.12 mas) from the 
finite source effect on the light curve, combined with an estimate of the angular size of the source given by the source 
position in a color-magnitude diagram. The lens mass, M — 0.612 ± 0.057 M , is found by combining these two 
quantities. This is the first time that parallax effects are detected for a caustic-crossing event and also the first time 
that the lens mass degeneracy has been completely broken through photometric monitoring alone. The combination 
of te and 8e also allows us to conclude that the lens lies in the near side of the disk, within 2.6 kpc of the Sun, while 
the radial velocity measurement indicates that the source is a Galactic bulge giant. 

Subject headings: gravitational lensing — binaries: general 



Although an initial objective of microlensing experiments was to probe the mass distribution of compact sub-luminous 
objects in the Galactic halo (Paczyhski 1986), the determination of the individual lens masses was in general believed not to 
be possible because the only generically measurable quantity is a degenerate combination of the mass with lens-source relative 
parallax and proper motion. However, Gould (1992) showed that there are additional observables that in principle can be used to 
break the degeneracy between physical parameters and so yield a measurement of the lens mass. If one independently measures 
the angular Einstein radius, 6*e, and the projected Einstein radius, te: 



'Department of Astronomy, the Ohio State University, 140 West 18th Avenue, Columbus, OH 43210, USA 

2 Department of Physics & Astronomy, University of Canterbury, Private Bag 4800, Christchurch, New Zealand 

3 Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD 21218, USA 

4 Institut d'Astrophysique de Paris, INSU-CNRS, 98 bis Boulevard Arago, F 75014 Paris, France 

5 South African Astronomical Observatory, P.O. Box 9, Observatory, 7935 South Africa 

6 Kapteyn Institute, Rijksuniversiteit Groningen, Postbus 800, 9700 AV Groningen, The Netherlands 

7 School of Natural Sciences, Institute for Advanced Study, Einstein Drive, Princeton, NJ 08540, USA 

8 Hubble Fellow 

'Physics Department, University of Tasmania, G.P.O. 252C, Hobart, Tasmania 7001, Australia 
'"School of Physics & Astronomy, University of St. Andrews, North Haugh, St. Andrews, Fife KY16 9SS, UK 
' ' Perth Observatory, Walnut Road, Bickley, Western Australia 6076, Australia 

12 Physics Department, Gettysburg College, 300 North Washington Street, Gettysburg, PA 17325, USA 



1. Introduction 
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then the lens mass M is obtained by 



M = 77^ E E 



= 0.1227 M, 



( 



1 AU 



) 




1 mas 



) 



(2) 



Here i? sc h = IGMc 2 is the Schwarzschild radius of the lens mass, D IC ] = D h 1 - D s 1 , and D L and D s are the distances to 
the lens and the source. 

Measurement of $e requires a standard ruler on the plane of the sky to be compared with the Einstein ring. For the subset of 
microlensing events in which the source passes very close to or directly over a caustic (region of singularity of the lens mapping 
at which the magnification for a point source is formally infinite), the finite source size affects the light curve. The source radius 
in units of the Einstein radius, p*, can then be determined through analysis of the light curve, and thus, one can determine 
(= 0*/p*) once the angular radius of the source, 9*, is estimated from its de-reddened color and apparent magnitude. Although 
this idea was originally proposed for point-mass lenses, which have point-like caustics (Gould 1994; Nemiroff & Wickramasinghe 
1994; Witt & Mao 1994), in practice it has been used mainly for binary lenses, which have line-like caustics and hence much 
larger cross sections (Alcock et al. 1997, 2000; Albrow et al. 1999a, 2000a, 2001a; Afonso et al. 2000). 

All ideas proposed to measure f E are based on the detection of parallax effects (Refsdal 1966; Grieger, Kayser, & Refsdal 
1986; Gould 1992, 1995; Gould, Miralda-Escude, & Bahcall 1994; Hardy & Walker 1995; Holz & Wald 1996; Honma 1999; An 
& Gould 2001): either measuring the difference in the event as observed simultaneously from two or more locations, or observing 
the source from a frame that accelerates substantially during the course of the event. The prime example of the latter is the Earth's 
orbital motion (annual parallax). For these cases, it is convenient to re-express r E as the microlens parallax, 7~e, 



Here, 7tl and ns are the (annual) trigonometric parallax of the lens and the source so that 7rx/l-mas = 1-kpc/Dx, where 
subscript 'X' is either 'L'(ens) or 'S'(ource). To date, 7te has been measured for several events by detecting the distortion in the 
light curve due to the annual parallax effect (Alcock et al. 1995; Mao 1999; Soszyhski et al. 2001; Bond et al. 2001; Smith, Mao, 
& Wozniak 2001; Mao et al. 2002; Bennett et al. 2001) 

An & Gould (2001) argue that 7r E is measurable for a caustic-crossing binary event exhibiting a well-observed peak caused 
by a cusp approach in addition to the two usual caustic crossings. Hence, one can determine the lens mass for such an event, 
since #e can be estimated from any well-observed caustic crossing. The caustic-crossing binary event EROS BLG-2000-5 is 
archetypal of such events. In addition, the event has a relatively long time scale (t E ~ 100 days), which is generally favorable 
for the measurement of 7Te- 

In fact, EROS BLG-2000-5 features many unique characteristics not only in the intrinsic nature of the event but also in the 
observations of it. These include a moderately well-covered first caustic crossing (entrance), a timely prediction - of not only the 
timing but also the duration - of the second caustic crossing (exit), the unprecedented four-day length of the second crossing, and 
two time-series of spectral observations of the source during the second crossing (Albrow et al. 2001b; Castro et al. 2001). To 
fully understand and utilize this wealth of information, however, requires detailed quantitative modeling of the event. Here, we 
present the first model of EROS BLG-2000-5 based on the photometric observations made by the Probing Lensing Anomalies 
NETwork (PLANET) 13 collaboration. In the current paper, we focus on the geometry of the event and find that both parallax 
and projected binary orbital motion are required to successfully model the light curve. Furthermore, we also measure the angular 
Einstein radius from the finite source effect during caustic crossings and the source angular size derived from the source position 
in a color-magnitude diagram (CMD). Hence, this is the first event for which both and #e are measured simultaneously and 
so for which the lens mass is unambiguously measured. 



7TE = 



7TL - 7TS _ l-kpc/A-d _ 1 AU 

#e 6*E/l-mas te 



(3) 
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2. Data 



On 2000 May 5, the Experience de Recherche d'Objets Sombres (EROS) collaboration issued an alert that EROS BLG- 
2000-5 was a probable microlensing event (RA = 17 h 53 m ll!5, Dec = -30°55'35"; I = 359?14, b = -2?43). On 2000 
June 8, the Microlensing Planet Search (MPS) 15 collaboration issued an anomaly alert, saying that the source had brightened by 
0.5 mag from the previous night and had continued to brighten by 0.1 mag in 40 minutes. PLANET intensified its observations 
immediately, and has been continuing to monitor the event up to the present (2001 August). Observations for PLANET were 
made from four telescopes: the Canopus 1 m near Hobart, Tasmania, Australia; the Perth/Lowell 0.6 m at Bickley, Western 
Australia, Australia; the Elizabeth 1 m at the South African Astronomical Observatory (SAAO) at Sutherland, South Africa; and 
the Yale/AURA/Lisbon/OSU (YALO) 1 m at the Cerro Tololo Interamerican Observatory (CTIO) at La Serena, Chile. Data were 
taken in V (except Perth), I (all four sites), and H (SAAO and YALO) bands. For the present study, we make use primarily of the 
/ band data. That is, we fit for the model using only / band data, while the V band data are used (together with the / data) only 
to determine the position of the source on CMD (see § 5). The light curve (Fig. 1) exhibits two peaks that have the characteristic 
forms of an entrance (A) and exit (B) caustic crossing (Schneider, Ehlers, & Falco 1992, also see fig. 1 of Gould & Andronov 
1999) immediately followed by a third peak (C) which is caused by the passage of the source close to a cusp. 

The data have been reduced in a usual way and the photometry on them has been performed by point spread function (PSF) 
fitting using DoPHOT (Schechter, Mateo, & Saha 1993). The relative photometric scaling between the different telescopes is 
determined as part of the fit to the light curve which includes the independent determinations of the source and the background 
fluxes at each telescope. We find, as was the case for several previous events, that due to the crowdedness of the field, the 
moderate seeing conditions, and possibly some other unidentified systematics, the amount of blended light entering the PSF is 
affected by seeing, and that the formal errors reported by DoPHOT tend to underestimate the actual photometric uncertainties. 
We tackle these problems by incorporating a linear seeing correction for the background flux and rescaling the error bars to force 
the reduced x 2 of our best model to unity. For details of these procedures, see Albrow et al. (2000a,b, 2001a) and Gaudi et al. 
(2002). 

We also analyze the data by difference imaging, mostly using ISIS (Alard 2000). We compare the scatter of the photometry 
on the difference images to that of the direct PSF fit photometry by deriving the normalized-summed squares of the signal-to-noise 
ratios: 



where F s is the source flux derived from the model, Ai is the magnification predicted by the model for the data point, Oi and <7j 
are the photometric uncertainty (in flux) for individual data before and after rescaling the error bars [i.e. di = cti(x 2 /N) 1 ! 2 ], 
and N is the number of the data points. Strictly speaking, Q defined as in equation (4) is model-dependent, but if the data sets 
to be compared do not differ with one another in systematic ways (and the chosen model is close enough to the real one), Q 
can be used as a proxy for the relative statistical weight given by the data set without notable biases. We find that difference 
imaging significantly improves the stability of the photometry for the data from Canopus and Perth, but it somewhat worsens 
the photometry for the data from SAAO and YALO. We suspect that the result is related to the overall seeing condition for the 
specific site, but a definite conclusion will require more detailed study and would be beyond the scope of the current paper. We 
hope to be able to further improve the photometry on difference images in the future, but for the current analysis, we choose to 
use the data set with the better Q so that only for Canopus and Perth data sets, we replace the result of the direct PSF photometry 
with the difference imaging analysis result. 

For the final analysis and the results reported here, we have used only a "high quality" subset of the data. Prior to any attempt 
to model the event, we first exclude various faulty frames and problematic photometry reported by the reduction/photometry 
software. In addition, data points exhibiting large (formal) errors and/or poor seeing compared to the rest of the data from the 
same site are eliminated prior to the analysis. In particular, the thresholds for the seeing cuts are chosen at the point where 
the behavior of "seeing-dependent background" becomes noticeably non-linear with the seeing variation. The criterion of the 
seeing and error cut for each data set is reported in Table 1 together with other photometric information. In conjunction with 
the proper determination of the error rescaling factors, we also remove isolated outliers as in Albrow et al. (2001a). Following 



http://www-dapnia.cea.fr/Spp/Experiences/EROS/alertes.html 
http://bustard.phys.nd.edu/MPS/ 




(4) 



-4- 



these steps, the "cleaned high-quality" /-band data set consists of 1286 (= 403 SAAO + 333 Canopus + 389 YALO + 161 Perth) 
measurements made during the 2000 season (between May 11 and November 12) plus 60 additional observations (= 25 SAAO 
+ 35 YALO) made in the 2001 season. Finally, we exclude 49 data points (= 19 SAAO + 20 Canopus + 10 Perth) that are very 
close to the cusp approach [2451736.8 < Heliocentric Julian Date (HJD) < 2451737.6] while we fit the light curve. We find 
that the limited numerical resolution of the source, which in turn is dictated by computational considerations, introduces errors 
in the evaluation of x 2 in this region of the order of a few, and in a way that does not smoothly depend on the parameters. These 
would prevent us from finding the true minimum, or properly evaluating the statistical errors. However, for the final model, we 
evaluate the predicted fluxes and residuals for these points. As we show in § 4, these resdiuals do not differ qualitatively from 
other residuals to the fit. 



3. Parameterization 

To specify the light curve of a static binary event with a rectilinear source trajectory requires seven geometric parameters: 
d, the binary separation in units of 9e", q, the binary mass ratio; a, the angle of the source-lens relative motion with respect to the 
binary axis; £e, the Einstein timescale (the time required for the source to transit the Einstein radius); uo, the minimum angular 
separation between the source and the binary center - either the geometric center or the center of mass - in units of #e; to, the time 
at this minimum; p*, the source size in units of #e- (In addition, limb-darkening parameters for each wave band of observations, 
and the source flux F s and background flux for each telescope and wavelength band are also required to transform the light 
curve to a specific photometric system.) Most generally, to incorporate the annual parallax and the projection of binary orbital 
motion into the model, one needs four additional parameters. However, their inclusion, especially of the parallax parameters, is 
not a trivial procedure, since the natural coordinate basis for the description of the parallax is the ecliptic system while the binary 
magnification pattern possesses its own preferred direction, i.e. the binary axis. In the following, we establish a consistent system 
to describe the complete set of the eleven geometric parameters. 



3.1. Description of Geometry 

First, we focus on the description of parallax. On the plane of the sky, the angular positions of the lens and the source (seen 
from the center of the Earth) are expressed generally by, 

0) = Vs,c + (t - *c) Ms + nss(t) , (5a) 

Vl(*) = ¥»l 1 c + (*-*c)Ml + ^(*)- (5b) 
Here ip S c and (p L c are the positions of the lens and the source at some reference time, t = t c , as they would be observed 
from the Sun, fi s and /x L are the (heliocentric) proper motion of the lens and the source, and s(t) is the Sun's position vector 
with respect to the Earth, projected onto the plane of the sky and normalized by an astronomical unit (see Appendix A). At any 
given time, t, s(t) is completely determined with respect to an ecliptic coordinate basis, once the event's (ecliptic) coordinates 
are known. For example, in the case of EROS BLG-2000-5 (A = 268?53, (3 = -7?50), <; = (0, r® sin/3) at approximately 2000 
June 19, where r® is the distance between the Sun and the Earth at this time in astronomical units (r® = 1.016). Then, the 
angular separation vector between the source and the lens in units of 8e becomes 

U(t) = ^ ~ Vh = V + (t- t c )n E - 7TE?(t) , (6) 

where v = (Vs.c^^Lcl/^E' Me = (Ms - MlV^e, and tte is defined as in equation (3). Although equation (6) is the most natural 
form of expression for the parallax-affected trajectory, it is convenient to re-express equation (6) as the sum of the (geocentric) 
rectilinear motion at the reference time and the parallactic deviations. In order to do this, we evaluate u and u at t = t c , 

u tc = u(t c ) =v- 7r E Q c , (7a) 

u tc = u(t c ) = Me - TEft, , (7b) 
where <; tc = <*(i c ) and <j tc = s(t c )- Solving equations (7) for /x E and v and substituting them into equation (6), one obtains 



u(t) = u tr +(t- t c )ii tc - 7r E £>p , 



(8) 
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where X>p = s(t) — q c — (t — t c )st c is the parallactic deviation. Note that T>p ~ (£t c /2)(i — i c ) 2 for t ~ i c , i.e. on relatively 
short time scales, the effect of the parallax is equivalent to a uniform acceleration of — 7TeQ c . Equation (8) is true in general for 
any microlensing event including point-source/point-lens events. 

Next, we introduce the binary lens system. Whereas the parallax-affected trajectory (eq. [8]) is most naturally described in 
the ecliptic coordinate system, the magnification pattern of the binary lens is specified with respect to the binary axis. Hence, 
to construct a light curve, one must transform the trajectory from ecliptic coordinates to the binary coordinates. If the origins of 
both coordinates are chosen to coincide at the binary center of mass, this transformation becomes purely rotational. Thus, this 
basically adds one parameter to the problem: the orientation of the binary axis in ecliptic coordinates. In accordance with the 
parameterization of the projected binary orbital motion, one may express this orientation using the binary separation vector, d, 
whose magnitude is d and whose direction is that of the binary axis (to be definite, pointing from the less massive to the more 
massive component). With this parameterization, the projection of the binary orbital motion around its center of mass is readily 
facilitated via the time variation of d. If the time scale of the event is relatively short compared to the orbital period of the binary, 
then rectilinear relative lens motion, d = d tc + d tc (t — t c ), will be an adequate representation of the actual variation for most 
applications (see e.g. Albrow et al. 2000a). Then, the light curve of a rotating binary event with parallax is completely specified 
by eleven independent parameters: d tc , d tc , u tc , iit c , q, p*, and 7Te- However, one generally chooses to make t c an independent 
parameter, such as the time when u tc ■ iit c = 0. In that case, the eleven parameters become d tc , d tc , u tc , Ut c , q, p*, tte, and t c . 

Although the parameterization described so far is physically motivated, and mathematically both complete and straightfor- 
ward, in practice it is somewhat cumbersome to implement into the actual fit. Therefore, we re-formulate the above parameter- 
ization for computational purposes. For the analysis of EROS BLG-2000-5, we first choose the reference time t c as the time of 
the closest approach of the source to the cusp, and rotate the coordinate system so that the whole geometry is expressed relative 
to the direction of d tc , i.e. the binary axis at time t c (see Fig. 2). We define the impact parameter for the cusp approach, u c 
(= \u tc — ttcuspl), and set u c > when the cusp is on the right-hand side of the moving source. Then, u tc is specified by 
t' E (= |« tc | _1 ), the instantaneous Einstein timescale at time t c , and by a', the orientation angle of u tc with respect to d tc . In 
addition, we express d tc in a polar-coordinate form and use the approximation that both the radial component, d, and tangential 
component, u>, are constant. Under this parameterization, d corresponds to the rate of expansion (d > 0) or contraction (d < 0) 
of the projected binary separation while u> is the angular velocity of the projected binary-axis rotation on the plane of the sky. Fi- 
nally, we define the microlens parallax vector, 7Te, whose magnitude is 7Te and whose direction is toward ecliptic west (decreasing 
ecliptic longitude). In the actual fit, 7r E ii and 7Te,_l, the two projections of 7Te along and normal to d tc , are used as independent 
parameters. Table 2 summarizes the transformation from the set of fit parameters (d tc ,q,a' ,u c ,t' E ,t c , p*, tte,|| > ^E,±,d, lo) to the 
set of the physical parameters (rf tc , d tc , u tc ,Ut c ,q, P*,^e, t c ). 



3.2. Terrestrial Baseline Parallax 

In general, the Earth's spin adds a tiny daily wobble of order ~ R®/fE (eq. [B3]), where R® is the Earth's radius, to the 
source's relative position seen from the center of the Earth as expressed in equation (6). Since R® = 4.26 x 10~ 5 AU, this 
effect is negligible except when the spatial gradient over the magnification map is very large, e.g. caustic crossings or extreme 
cusp approaches. Even for those cases, only the instantaneous offsets are usually what matters because the source crosses over 
the region of extreme gradient with a time scale typically smaller than a day. Hence, unless the coverage of the crossing from 
two widely separated observers significantly overlaps, the effect has been in general ignored when one models microlensing light 
curves. 

However, in case of EROS BLG-2000-5, the second caustic crossing lasted four days, and therefore daily modulations 
of magnifications due to the Earth's rotation, offset according to the geographic position of each observatory, may become 
important, depending on the actual magnitude of the effect (See also Honma 1999 for a similar discussion on the short time 
scale magnification modulation observed from an Earth-orbiting satellite). Hardy & Walker (1995) and Gould & Andronov 
(1999) investigated effects of the terrestrial baseline parallax, for fold-caustic crossing microlensing events mainly focused on the 
instantaneous offsets due to the separation between observers. They argued that the timing difference of the trailing limb crossing 
for observations made from two different continents could be of the order of tens of seconds to a minute (Hardy & Walker 1995) 
and the magnifications near the end of exit-type caustic crossings could differ by as much as a few percent (Gould & Andronov 
1999). Suppose that (j> 2 is the angle at which the source crosses the caustic, A max is the magnification at the peak of the crossing, 
and A cc is the magnification right after the end of the crossing. Then, for the second caustic crossing of EROS BLG-2000-5, 
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since csc ^(Re/^E) — 9 (0.25/ sin 02 ) min and p~ 1 (A maK /A cc )(R^/f-E) — 2 x 10~ 2 , the time for the end of the second 
caustic crossing may differ from one observatory to another by as much as ten minutes and the magnification difference between 
them at the end of the crossing can be larger than one percent, depending on the relative orientation of observatories with respect 
to the event at the time of the observations. Based on a model of the event, we calculate the effect and find that it causes the 
magnification modulation of an amplitude as large as one percent (Fig. 3). In particular, night portions of the observatory-specific 
light curves exhibit steeper falls of flux than would be the case if the event were observed from the Earth's center. That is, the 
source appears to move faster during the night because the reflex of the Earth's rotation is added to the source motion. This would 
induce a systematic bias in parameter measurements if it were not taken into account in the modeling. We thus include the (daily) 
terrestrial baseline parallax in our model to reproduce the observed light curve of EROS BLG-2000-5. Here, we emphasize that 
this inclusion requires no new free parameter for the fit once geographic coordinates of the observatory is specified and R® in 
units of AU is assumed to be known (see Appendix B). 

In addition, we note the possibility of a simple test of the terrestrial baseline parallax. Figure 4 shows that the end of the 
crossing observed from S AAO is supposed to be earlier than in the geocentric model - and earlier than seen from South American 
observatories. Unfortunately, near the end of the second caustic crossing (HJD ~ 2451733.66), PLANET data were obtained 
only from SAAO near the very end of the night - the actual trailing limb crossing is likely to have occurred right after the end of 
the night at SAAO, while YALO was clouded out due to the bad weather at CTIO (The event was inaccessible from telescopes 
on Australian sites at the time of trailing limb crossing.) However, it is still possible to compare the exact timing for the end of the 
second crossing derived by other observations from South American sites with our model prediction and/or the observation from 
SAAO. In particular, the EROS collaboration has published a subset of their observations from the Marly 1 m at the European 
Southern Observatory (ESO) at La Silla, Chile, for the second crossing of EROS BLG-2000-5 (Afonso et al. 2001). Comparison 
between their data and our model/observations may serve as a confirmation of terrestrial parallax effects. 



3.3. Limb-darkening Coefficients 



Because of the unprecedentedly long time scale of the second caustic crossing and the extremely close approach to the cusp, 
as well as the high quality of the data, we adopt a two-parameter limb-darkening law of the form, 



(1 - T A - A A ) + ^ cos + ^ cos 1 / 2 # 



(9) 



to model the surface brightness profile of the source Here S\ = F Si \/(w6l) is the mean surface brightness of the source and d is 
the angle between the normal to the stellar surface and the line of sight, i.e., siwd = 6/0* where 9 is the angular distance to the 
center of the source. This is an alternative form of the widely-used square-root limb-darkening law, 



S\{&) = S x (0) \l - c A (l - cos i?) - d A (l - Vcosd) 



(10) 



However, instead of being normalized to have the same central intensity S A (0) as a uniform source, the form we adopt (eq. [9]) 
is normalized to have the same total flux F s A = (27r# 2 ) 5 A (i?) sini? d(sin#). That is, there is no net flux associated with the 
limb-darkening coefficients. The transformation of the coefficients in equation (9) to the usual coefficients used in equation (10) 
is given by 

6r A , 5A A 



CA = 



dx = 



(11) 



4 + 2T A + A A ' " 4 + 2T A +A A 
Note that, although we fit limb-darkening, the detailed discussion of its measurement and error analysis will be given elsewhere. 



4. Measurement of the Projected Einstein Radius 

Table 3 gives the parameters describing the best-fit microlens model (see Appendix C for details of modeling and Ap- 
pendix D for the discussion on the error determination) for the PLANET /-band observations of EROS BLG-2000-5. We also 
transform the fit parameters to the set of parameters introduced in § 3.1. In Table 4, we provide the result of limb-darkening 
coefficient measurements. The measurements of two coefficients, Tj and Aj are highly anti-correlated so that the uncertainty 
along the major axis of error ellipse is more than 50 times larger than that along its minor axis. While this implies that the 
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constraint on the surface brightness profile derived by the caustic crossing light curve is essentially one dimensional, its natural 
basis is neither the linear nor the square-root parameterized form. We plot, in Figure 5, the surface brightness profile indicated 
by the fit and compare these with theoretical calculations taken from Claret (2000). The figure shows that allowing the profile 
parameters to vary by 2-cr does not have much effect on the central slope, but cause a large change in the behavior near the limb. 
We speculate that this may be related to the specific form of the time sampling over the stellar disk, but more detailed analysis 
regarding limb darkening is beyond the scope of the present paper and will be discussed elsewhere. 

Figures 6 and 7 show "magnitude residuals", (2.5/ In 10)[AF / (AFs)], from our best model. Note in particular that the 
residuals for the points near the cusp approach (HJD = 2451736.944) that were excluded from the fit do not differ qualitatively 
from other residuals. It is true that the mean residual for Perth on this night (beginning 4.0 hrs after t c and lasting 2.3 hrs) was 
about 2% high. However, the Canopus data, which span the whole cusp approach from 1.2 hrs before until 7.7 hrs after t c , agree 
with the model to within 0.5%. Moreover, the neighboring SAAO and YALO points also show excellent agreement. See also 
Fig. 8 and especially Fig. 9. 

From the measured microlens parallax (tte — 0.277 ± 0.008), the projected Einstein radius is (eq. [3]); 

r E = = 3.61 ±0.11 AU . (12) 

We also derive £e, the heliocentric Einstein timescale; 

Me = «t c + ^E^tc , (13a) 
t E = I/XeT 1 = 76.8 ±2.1 days, (13b) 

and the direction of /x E is 36? 3 from ecliptic west to south. Since, toward the direction of the event, the Galactic disk runs along 
60?2 from ecliptic west to south, fi E points to 23?9 north of the Galactic plane. The overall geometry of the model is shown in 
Figure 8, and Figure 9 shows a close-up of the geometry in the vicinity of the cusp approach (t ~ t c ). Next, the projected velocity 
v on the observer plane is 

v = r EfJ , E = (-74.5 ± 3.1, 33.0 ±11.0) km s" 1 . (14) 

Here the vector is expressed in Galactic coordinates. We note that the positive x-direction is the direction of the Galactic rotation 
- i.e. the apparent motion of the Galactic center seen from the Local Standard of Rest (LSR) is in the negative x direction - while 
the positive y-direction is toward the north Galactic pole so that the coordinate basis is left-handed. 



5. Measurement of the Angular Einstein Radius 

The angular radius of the source, 6>», is determined by placing the source on an instrumental CMD and finding its offset 
relative to the center of the red giant clump, whose de-reddened and calibrated position toward the Galactic bulge are known from 
the literature. For this purpose, we use the data obtained from YALO. We find the instrumental I magnitude of the (de-blended) 
source I s from the fit to the light curve. To determine the color, we first note that, except when the source is near (and so partially 
resolved by) a caustic, microlensing is achromatic. That is, the source is equally magnified in V and /: Fy = Fy.b + AF v, s and 
Fi = Fj ,b + AFi tS . Hence, we assemble pairs of V and / points observed within 30 minutes (and excluding all caustic-crossing 
and cusp-approach data) and fit these to the form Fy = a\ + a 2 Fj. We then find (V — /) s ,inst = —2.5 loga 2 . We also find the 
magnitude of the blend Jb from the overall fit, and solve for the color of the blend (V — /)b,inst = —2.5 log[a 2 + (a\/Fj^)]. 

We then locate the source on the instrumental CMD and measure its offset from the center of the red giant clump; A(V—I) = 
0.276 ± 0.010 and AI = 0.33 ± 0.02. Here, the quoted uncertainty includes terms for both the clump center and the source 
position. Finally, using the calibrated and de-reddened position of the red clump, [(V — I)o, Io] = [1-114± 0.003, 14.37± 0.02] 
(Paczyriski et al. 1999), the source is at (V - J) = 1.390 ± 0.010, J = 14.70 ± 0.03. The procedure does not assume any 
specific reddening law for determining the de-reddened source color and magnitude, but may be subject to a systematic error due 
to differential reddening across the field. 

We also perform a photometric calibration for observations made at SAAO. The calibration was derived from images ob- 
tained on three different nights by observing Landolt (1992) standards in the 2000 season, and the resulting transformation equa- 
tion were confirmed with a further night's observations in 2001. Applying the photometric transformation on the fit result for the 
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source flux, we obtain the - unmagnified - source magnitude of I = 16.602 in the standard (Johnson/Cousins) system, while the 
standard source color is measured to be V — I = 2.69. By adopting the photometric offset between SAAO and YALO instrumental 
system determined by the fit, we are able to derive a "calibrated" YALO CMD for the field around EROS BLG-2000-5 (Fig. 10). 
The positions of the source (S), blend (B), and clump center (C) are also overlaid on the CMD. The CMD also implies that the 
reddening toward the field is A/ = 1.90 and E(V - I) = 1.30, which yields a reddening law, R VI = A v / ' E(V - I ) = 2.46, 
nominally consistent with the generally accepted Ryi = 2.49 ± 0.02 (Stanek 1996). If we adopt a power-law extinction model 
A x oc \- p , then R VI = 2.46 corresponds to p = 1.37. For this index, E(B-V)/E(V - I) = 0.838, which predicts a lower 
extinction than the spectroscopically determined reddening of E(B — V) = 1.30 ± 0.05 by Albrow et al. (2001b). 

Both the source and the blend lie redward of the main population of stars in the CMD. One possible explanation is that the 
line of sight to the source is more heavily reddened than average due to differential reddening across the field. Inspection of the 
images does indicate significant differential reddening, although from the CMD itself it is clear that only a small minority of stars 
in the field could suffer the additional AA V ^0.7 that would be necessary to bring the source to the center of the clump. There is 
yet another indication that the source has average extinction for the field - i.e. the same or similar extinction as the clump center. 
The de-reddened (intrinsic) color derived on this assumption, (V — I)q = 1.39 is typical of a K3III star (Bessell & Brett 1988), 
which is in good agreement with the spectral type determined by Albrow et al. (2001b). On the other hand, if the source were a 
reddened clump giant with (V — 7)o ~ 1.11, then it should be a Kl or K2 star. 

We then apply the procedure of Albrow et al. (2000a) to derive the angular source radius 0* from its de-reddened color and 
magnitude: first we use the color-color relations of Bessell & Brett (1988) to convert from (V — I)o to (V — K)q, and then we 
use the empirical relation between color and surface brightness to obtain 0* (van Belle 1999). From this, we find that 

9* = 6.62 ± 0.58 ^as 

= 1.42 ± 0.12 Rokpc- 1 , (15) 

where the error is dominated by the 8.7% intrinsic scatter in the relation of van Belle (1999). Alternatively, we use the different 
calibration derived for K giants by Beuermann, Baraffe, & Hauschildt (1999) and obtain 6* = 6.47 /zas, which is consistent 
with equation (15). Finally, we note that if the source were actually a clump giant that suffered AAy = 0.7 extra extinction, 
its angular size would be 7.08 ^as. Since we consider this scenario unlikely, and since in any event the shift is smaller than the 
statistical error, we ignore this possibility. 

From this determination of 0* and the value, p* = (4.80 ± 0.04) x 10~ 3 , determined from the fit to the light curve, we 
finally obtain, 

6 E = — = 1.38 ±0.12 mas, (16) 

P* 

where the error is again dominated by the intrinsic scatter in the relation of van Belle (1999). 



6. The Lens Mass and the Lens-Source Relative Proper Motion 

By combining the results obtained in §§ 4 and 5, one can derive several key physical parameters, including the lens mass 
and the lens-source relative parallax and proper motion; 

c 2 

M = — f E 6» E = 0.612 ±0.057 M , (17) 

7J rel = -f- = 2.62 ±0.24 kpc; (18a) 
7Tl — 7Ts = = 0.382 ± 0.035 mas . (18b) 

A-el 

I Ms ~ Ml I = 0e|MeI = 18.0 ± 1.7 ^as day -1 

= 31.1 ±2.9 km s- 1 kpc" 1 . (19) 
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For the binary mass ratio of the best-fit model, q = 0.749, the masses of the individual components are 0.350 M and 
0.262 M Q , both of which are consistent with the mass of typical mid-M dwarfs in the Galactic disk. Using the mass-luminosity 
relation of Henry & McCarthy (1993), and adopting M v = 2.89 + 3.37(V - I), the binary then has a combined color and 
magnitude Mj = 8.2 ± 0.2 and (V - I) = 2.54 ± 0.08. Since D~^ = 1 - D s \ £> rel is an upper limit for D L , i.e. the lens is 
located in the Galactic disk within 2.6 kpc of the Sun. Furthermore, an argument based on the kinematics (see § 7) suggests that 
Dg > 8 kpc so that it is reasonable to conclude that <~ 2 kpc (m — M ~ 11.5). Hence, even if the binary lens lay in front of 
all the extinction along the line of sight, it would be ~ 2 mag fainter than the blend (B), and so cannot contribute significantly to 
the blended light. However, the proximity of the lens to the observer is responsible for the event's long time scale and quite large 
parallax effect. 



7. The Kinematic Constraints on the Source Distance 

The projected velocity v (eq. [14]) is related to the kinematic parameters of the event by 

* = D rcl (n s - fijj 



= (# re] g - D rel g- + v Q J , (20a) 

where the final subscript, _L, serves to remind us that all velocities are projected onto the plane of the sky. Writing v Q — 
Viot + v q,p and = v Iot + t>L, P as the sum of the Galactic rotation and the peculiar velocities and eliminating £>l in favor of 
-D re i and Dg, equation (20a) can be expressed as 



v = 



D Te \, ( D re l . 

{V S ~ «rot) - 1 + -=— t>L,p + f0,p 



D S V D S 



(20b) 



For a fixed Ds and with a known Galactic kinematic model, the distribution function of the expected value for v can be derived 
from equation (20b). The measured value of v can therefore be translated into the relative likelihood, C for a given Ds and the 
assumed kinematic model, 

-21n£ = (v- (v)) ■ C" 1 • (v - (v)) + In |C| + constant, (21a) 

(v) = X S ((v S ,±) - Urot) - (1 + Xs)(v L . p .±) + U0,p,± , (21b) 

C = x 2 s C[v s ,±] + (1 + x s ) 2 C[v L ,±] + C[v] , (21c) 

where C[t?s,x] an d are the dispersion tensors for the source and the lens transverse velocity, C[v] is the covariance 

tensor for the measurement of v, and xs = D rc \/Ds- 

We evaluate £ as a function of Ds assuming three different distributions for vs.±_ which respectively correspond to source 
being in the near disk ((v s .±) = v Iot + («s, P ,J.». bulge ((v s .±) = («s, P ,J.». and the far disk ((v s .±) = -v mt ~ (v s , p ,±)). 
Adopting the Galactic rotational velocity, u rot = (220, 0) km s _1 , the solar motion, i>©,p,_L = (5.2, 7.2) km s _1 (Binney & 
Merrifield 1998), and the kinematic characteristics of the lens and the source given in Table 5, we derive the relative likelihood 
as a function of source position (Fig. 11), and find that the measured value of v (eq. [14]) mildly favors the far disk over the 
bulge as the location of the source by a factor of <~ 2.3. The near disk location is quite strongly disfavored (by a factor of <~ 10.6 
compared to the far disk, and by a factor of ~ 4.7 compared to the bulge). 

The complete representation of the likelihood for the source location requires one to consider all the available constraints 
relevant to the source distance. In particular, these include the radial (line-of-sight) velocity measurement and the number density 
of stars constrained by the measured color and brightness (or luminosity) along the line of sight. Although one may naively 
expect that a bulge location is favored by the high density of stars in the bulge, which follows from the fact that the line of sight 
passes within 400 pc of the Galactic center, it is not immediately obvious how the disk contribution compares to the bulge stars 
for the specific location of the source on the CMD (Fig. 10). Even with no significant differential reddening over the field, the 
particular source position, which is both fainter and redder than the center of the clump, can be occupied by relatively metal-rich 
(compared to the bulge average) first ascent giants in either the bulge or the far disk. Because first ascent giants become redder 
with increasing luminosity, the source must have higher metallicity if it lies in the bulge than in the more distant disk. There 
exist two estimate of the source metal abundance: [Fe/H]= —0.3 ± 0.1 by Albrow et al. (2001b) based on VLT FORS1 spectra 
and [Fe/H]= —0.3 ± 0.3 by Minniti et al. (2002) based on VLT UVES spectra. However, to incorporate these measurement into 
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a likelihood estimate would require a more through understanding of all the sources of error as well as a detailed model of the 
metallicity distribution of the bulge and the disk. For this reason, we defer the proper calculation of the likelihood in the absence 
of any definitive way to incorporate the specific density function, and restrict ourselves to the kinematic likelihood. 

We determine the radial velocity of the source from Keck HIRES spectra of Castro et al. (2001), kindly provided to us by R. 
M. Rich. We find the line-of-sight velocity to be <~ —100 km s _1 (blueshift; heliocentric) 16 , which translates to ~ —90 km s _1 
Galactocentric radial velocity accounting for the solar motion of 10 km s _1 . The derived radial velocity strongly favors bulge 
membership of the source since it is three times larger than the line-of-sight velocity dispersion for disk K3 giants, but is consistent 
with the motions of typical bulge stars (c.f. Table 5). Because the correlation between the radial and the transverse velocity for 
K3 giants either in the disk or in the bulge is very small, the likelihood for the radial velocity can be estimated independently 
from the likelihood for the transverse velocity, and the final kinematic likelihood is a simple product of two components. We 
find that the final kinematic likelihood including the radial velocity information indicates that the source star belongs to the bulge 
(Fig. 1 1): the likelihood for the bulge membership is about six times larger than that for the far disk membership. 

Finally, we note that, if the source lay in the far disk, it should have an additional retrograde proper motion of <~ 4 mas yr~ 1 x 
(10 kpc/Ds) with respect to the bulge stars, which should be measurable using adaptive optics or the Hubble Space Telescope. 



8. Consistency between the Lens Mass and the Binary Orbital Motion 

For the derived parameters, we find a projected binary lens separation r± = dD^O^ = 5.52 AU (5.25 AU) and transverse 
orbital speed v ± = [d 2 + (uod) 2 ] 1 / 2 D L d E = 2.76 km s" 1 (2.62 km s" 1 ) for D s = 10 kpc (8 kpc). We can now define the 
transverse kinetic and potential energies Kj_ = [<?/(! + q) 2 ]Mv\/2 and Tj_ = —[q/(l + q) 2 ]GM 2 /r±, and evaluate their ratio, 



K± 



2GM c 2 rM^ 1 + (D S 8 E )-^ 



i: - ' 2 5.9 : 5.1 (D s 10 kpc! . (22) 



r±v 2 ± 2 d[d, 2 + {uod) 2 ] 

Here, the error is dominated by the uncertainty in the measurement of ui. Since r± < r and v± < v, it follows that \T± | > |T| and 
K± < K, and hence \T±/K±\ > \T/K\, where \T/K\ is the ratio of the true three-dimensional potential and kinetic energies. 
The latter must be greater than unity for a gravitationally bound binary. This constraint is certainly satisfied by EROS BLG-2000- 
5. Indeed, perhaps it is satisfied too well. That is, what is the probability of detecting such a large ratio of projected energies? To 
address this question, we numerically integrate over binary orbital parameters (viewing angles, the orbital orientation and phase, 
and the semi-major axis) subject to the constraint that r± = 5.5 AU (Ds — 10 kpc), and at several fixed values of the eccentricity, 
e. We assume a random ensemble of viewing angles and orbital phases. The results shown in Figure 12 assume a Duquennoy & 
Mayor (1991) period distribution, but are almost exactly the same if we adopt a flat period distribution. All of the eccentricities 
shown in Figure 12 are reasonably consistent with the observed ratio, although higher eccentricities are favored. 

One can also show that high eccentricities are favored using another line of argument. First, note that \r j_ x v± | = lot\ is the 
same as the projection of the specific angular momentum of the binary to the line of sight, |(r x v) ■ n\ = 2ira 2 (l — e 2 ) \ cosi\/P. 
Here a and P are the semi-major axis and the orbital period of the binary, n is the line-of-sight unit vector, and i is the inclination 
angle of the binary. Combining this result with Kepler's Third Law, 47r 2 a 3 = GMP 2 , one obtains (here and throughout this 
section, we assume that Ds = 10 kpc); 

(cur 2 ) 2 4 U! 2 d 4 

a(l - e 2 ) 2 cos 2 i = y J- I = — = — — = 0.0013 ± 0.0341 AU . (23) 

GM c 2 r E E [fi + We)" 1 ] 4 

From the energy conservation, (v 2 /2) — (GM/ r) = —GM/ (2a), one may derive a lower limit for a; 

r ( rv 2 \ 1 r 
a =2{ 1 -2GM) > ^^1)= 2-87 ± 0.09 AU. (24) 

The corresponding lower limit for the binary period is P > 6.22 yr for Ds — 10 kpc. Now, we can derive a constraint on the 
allowed eccentricity and inclination by dividing equation (23) by equation (24); 

(l-e 2 )|cosi| < 0.021 ± 0.281. (25) 



16 Recently Minniti et al. (2002) reported a radial velocity measurement of —191 km s 1 for the source star of EROS BLG-2000-5. At this point, we do not 
know what is the reason for this discrepancy. 
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The constraint here essentially arises from the fact that the fit barely detects projected angular motion u>, while the formal precision 
of its measurement corresponds to <~ 80 yr orbital period in l-cr level. For the observed projected separation and the derived 
binary mass, this apparent lack of the binary angular motion therefore naturally leads us to conclude that the binary orbit is either 
very close to edge-on or highly eccentric (or both). 



9. Another Look at MACHO 97-BLG-41 

We have found that both parallax and binary orbital motion are required to explain the deviations from rectilinear motion 
exhibited by the light curve of EROS BLG-2000-5. In a previous paper about another event, MACHO 97-BLG-41 (Albrow et 
al. 2000a), we had ascribed all deviations from rectilinear motion to a single cause: projected binary orbital motion. Could both 
effects have also been significant in that event as well? Only detailed modeling can give a full answer to this question. However, 
we can give a rough estimate of the size of the projected Einstein ring that would be required to explain the departures from 
linear motion seen in that event. 

Basically what we found in the case of MACHO 97-BLG-41 was that the light curve in the neighborhood of the central 
caustic (near HJD ~ 2450654) fixed the lens geometry at that time, and so predicted both the position of the outlying caustics 
and the instantaneous velocity of the source relative to the Einstein ring. If this instantaneous relative velocity were maintained, 
then the source would have missed this outlying caustic by about Au obs ~ 0.4 (Albrow et al. 2000a, fig. 3). On the other hand, 
the predicted displacement of a caustic due to parallax is 

Ait prod = -7T E X>p ; (26a) 

T> P = St cc l - s toc , 2 - (i C c,i - *cc,2Kt cc , 2 • (26b) 

We find |X> P | = 0.072, and hence, 

7Te = 5.6 • (27) 

AUobs 

That is, to explain by parallax the order of the effect seen requires r-g ~ 0.18 AU, which (using the measured 9e = 0.7 mas) 
would in turn imply a lens mass M ~ 0.015 M Q , a lens distance of Z?l ^ 250 pc, and the projected lens-source relative transverse 
velocity on the observer plane (at time t cc ,2) of only 13 km s _1 . While these values cannot be strictly ruled out, they are a priori 
extremely unlikely because the optical depth to such nearby, low-mass, and slow lenses is extremely small. On the other hand, if 
the lens lies at a distance typical of bulge lenses ~ 6 kpc and the source is a bulge star, then ?e — -D to i$e ~ 17 AU, which 
would imply Aw pre d/Aw obs <~ 1%. That is, parallax would contribute negligibly to the observed effect. We therefore conclude 
that most likely parallax does not play a major role in the interpretation of MACHO 97-BLG-41, but that detailed modeling will 
be required to determine to what extent such a role is possible. 

What is the physical reason that parallax must be so much stronger (i.e., 7te must be so much bigger) to have a significant 
effect in the case of MACHO 97-BLG-41 than for EROS BLG-2000-5? Fundamentally, the former is a close binary, and there is 
consequently a huge "lever arm" between the radial position of the outlying caustic, u c <~ 1.7 and the binary separation, d <~ 0.5. 
That is u c /d ~ 3.2. This is almost an order of magnitude larger than for EROS BLG-2000-5, for which u c ~ 0.8 (Fig. 8), 
d <~ 1.9, and u c /d ~ 0.4. Consequently, parallax has to be an order of magnitude larger to reproduce the effects of the same 
amount of the projected binary orbital motion. 



10. Summary 

We have presented here results from two seasons of photometric /-band monitoring of the microlensing event EROS BLG- 
2000-5, made by the PLANET collaboration. The light curve exhibits two peaks which are well explained by a finite source 
crossing over a fold-type - inverse-square-root singularity - caustic, followed by a third peak which is due to the source's passage 
close to a cusp. We find no geometry involving a rectilinear source-lens relative trajectory and a static lens that is consistent 
with the photometry. However, by incorporating both parallax and binary orbital motion, we are able to model the observed 
light curve. In particular, the detection of the parallax effect is important because it enables us to derive the microlens mass, 
M = 0.612 ± 0.057 Mq, unambiguously by the combination of the projected Einstein radius - measured from the parallax effect 
- and the angular Einstein radius - inferred from the source angular size and the finite source effect on the light curve. The source 
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size is measurable from the magnitude and the color of the source. The kinematic properties of the lens/source system derived 
from our model together with the lens-source relative parallax measurement as well as an additional radial velocity measurement 
indicate that the event is most likely caused by a (K3) giant star in the bulge being lensed by a disk binary (M dwarf) system 
about 2 kpc from the Sun. Additional information on the specific density function along the line of sights, differential reddening 
across the field, and a metallicity measurement of the source could further constrain the source location. 
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A. Determination of ^ 

If s is the Sun's position vector with respect to the Earth normalized by 1 AU, then the projection of s onto the plane of the 

sky, is 

S = s — (s ■ h)n , (Al) 

where h is the line-of-sight unit vector toward the position of the event on the sky, while the projection of p, the unit vector 
toward the north ecliptic pole (NEP), is given by p = p — (p ■ n)h. Then, (<; w , <; n ), the ecliptic coordinate components of are 

_ (p x <r) ■ n _ (pxs)-n 

^ = £± = -<*;* ) ^%, ( A2b) 

\p\ y/l-(p- n) 2 

where we make use of p ■ s = 0. One can choose three-dimensional coordinate axes so that the x-axis is the direction of the 
vernal equinox, the z-axis is the direction to the NEP, and y = z x x. Then, 

p= (0,0,1); (A3a) 

s = (r e cosA ,r e sinA ,O) ; (A3b) 

n = (cos Ao cos /3o , sin Ao cos 0o , sin /3 ) , ( A3c) 

where r ffi is the distance to the Sun from the Earth in units of AU, A Q is the Sun's ecliptic longitude, and (Ao, /?o) are the ecliptic 
coordinates of the event. By substituting equations (A3) into equations (A2), one finds that 

Sw = -r® sin(A Q - A ) , (A4a) 

C„ = -r® cos(A Q - A ) s'mPo . (A4b) 



In general, one must consider the Earth's orbital eccentricity (e = 1.67 x 10 2 ) to calculate <; for any given time. Then, 

r ffi = 1 - e cos tp ; A Q = £ - <f>r , (A5a) 
ip — e sin ip = Clt , (A5b) 
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. > (1 — e 2 ) 1 ' 2 sin ib _ cosib — e 

sin£=^- ^ -r; cos£=- r, (A5c) 

1 — e cos ^ 1 — e cos tp 

where ^> and £ are the eccentric and true anomalies of the Earth, T = 77? 86 is the true anomaly at the vernal equinox (March 
20, 07 h 35 m UT for 2000; Larsen & Holdaway 1999), t is the time elapsed since perihelion, and Q = 2n yr _1 . Note that the 
Earth was at perihelion at January 3, 05 h UT for 2000 (Larsen & Holdaway 1999). Although equations (A5) cannot be solved 
for r ffi and A© in closed form as functions of t, one can expand in series with respect to e and approximate up to the first order 
(epicycle approximation) so that 

r© = 1 - e cos(Oi) ; A Q = fit - <fn + 2e sin(fii) . (A6) 



B. Microlens Diurnal Parallax 

The angular position of a celestial object observed from an observatory on the surface of the Earth is related to its geocentric 
angular position, y x by 

-. \{ 

¥>x = ¥>x ~ -J^- . (Bl) 

where 

7 = r — (r ■ h)n (B2) 

is the projection of the position vector, r, of the observatory with respect to the center of the Earth, onto the plane of the sky 
and normalized by the mean radius of the Earth, R ffi ; and Dx and h are the distance and the line-of-sight vector to the object 
of interest. If one observe a microlensing event, the actual dimensionless lens-source separation vector therefore differs from u 
(eq. [6]) by 

u> = ^-^ = M+ ^7 

= v + (t-t c )v E -M<;--lj^j) ■ (B3) 

To find the algebraic expression for ecliptic coordinate components of 7, we choose the same coordinate axes as for equa- 
tions (A3). Then, the position vector, r is expressed as 

r = (cos 6 g cos 7~x , cos 5 g sin rx cos e + sin S g sin e, sin S g cos e — cos S g sin rx sin e) , (B4) 

where <5 g is the geographic latitude of the observatory, rx is the hour angle of the vernal equinox (i.e. the angle of the local 
sidereal time) at the observation, and e = 23? 44 is the angle between the direction toward the north celestial pole and the NEP 
(here we also assume that the Earth is a perfect sphere). Then, following a similar procedure as in Appendix A, 

(px7) 1 n = iP-r)-n 
\p\ v 71 - (p-n) 2 

V ■ 1 P r (p-n)(r-n) 

7™ = pr"j = ===== — , (B5b) 

\P\ 0-^ - n ) 

one obtains the ecliptic coordinate components of 7, 

"fw = — sin S g sin e cos Ao + cos 5 g (cos rx sin Ao — sin rx cos e cos A ) ; (B6a) 
7„ = sin Sg (cose cos (3q — sine sin Aq sin/3 ) — cos<5 g [cosrx cos Aq sin/3 + sin rx (sin e cos (3q + cose sin Aq sin/3 )] ■ (B6b) 



C. The Choice of Fit Parameters 

Judged by the number of fitting parameters alone, EROS BLG-2000-5 is by far the most complex event ever analyzed: 
compared to the runner-up, MACHO 97-BLG-41 (Bennett et al. 1999; Albrow et al. 2000a), it has two more geometric parameters 
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and one more limb-darkening parameter. As a direct result, the path toward choosing a modeling procedure was substantially 
more tortuous than usual. We therefore believe that it is important to document this path, at least in outline, in order to aid in the 
modeling of future events. 

As stated in § 3, the seven standard parameters for binary events are (d, q, a, u , tg, to, P*)- Immediately following the 
first caustic crossing, we fit this crossing to five empirical parameters, including t CCi i and Ati, the time and half-duration of the 
first crossing. We then changed our choice of parameters (to, p») — > (t CCj i, Ati) according to the prescription of Albrow et al. 
(1999b), effectively cutting the search space down from seven to five dimensions, and speeding up the search accordingly. This 
permitted us to accurately predict in real time not only the time, but also the (4-day) duration of the second crossing which in turn 
allowed two groups to obtain large-telescope spectra of the crossing (Albrow et al. 2001b; Castro et al. 2001). This was the first 
prediction of the duration of a caustic crossing. 

Why is the substitution (to, p*) — ► (t C c,i, Ati) critical? Both t CCi i and Ati are determined from the data with a precision 
~ 10~ 3 day ~ 10~ 5 t^. Hence, if any of the parameters, a, u , to, p*, are changed individually by 10~ 4 (subsequent changes 
of t CC; i and Ati by 10~ 4 t^), this will lead to an increase A\ 2 ~ 100. As a result, even very modest movements in parameter 
space must be carefully choreographed to find a downhill direction on the \ 2 hypersurface. By making (t cc ,i, Ati) two of the 
parameters and constraining them to very small steps consistent with their statistical errors, one in effect automatically enforces 
this choreography. 

In all the work reported here, we searched for \ 2 minima at fixed (d, q), and repeated this procedure over a (d, q) grid. We 
found for this event (as we have found for others) that regardless of what minimization technique we apply, if we search (d, q) 
space simultaneously with the other parameters, then either we do not find the true \ 2 minimum or the search requires prohibitive 
amounts of time. 

Following the second crossing, we added a linear limb-darkening (LD) parameter, but otherwise continued with the same 
parameterization. We found that the fitting process was then enormously slowed down because small changes a, uo, or t^ led to 
large changes in t CCt2 and At 2 (the second crossing time and half-duration), whereas these quantities were directly fixed by the 
data. We therefore changed parameters (a, uo, t-&) — > (t C c,2, At 2 , t ax is), where t ax ; s is the time the source crossed the cusp axis. 
Hence, all five of the non-(d, q) parameters were fixed more-or-less directly by the data, which greatly improved the speed of our 
parameter search. We thus quickly found the x 2 minimum for this (7+l)-parameter- seven geometry plus one limb-darkening - 
space. (Note that the LD parameter, like the source flux, the background flux, and the seeing correlation, is determined by linear 
fitting after each set of other parameters is chosen. Hence, it exacts essentially zero computational cost. We therefore track it 
separately.) 

Since the light curve showed systematic residuals of several percent (compared to daily-averaged photometry errors <C 1%), 
we were compelled to introduce more parameters. We first added two parallax parameters; the magnitude 7Te and its relative 
orientation with respect to the binary, yielding a (9+l)-parameter fit. Since five of the nine geometrical parameters remained 
empirical, this procedure also converged quickly. However, while \ 2 had fallen by several hundred (indicating a very significant 
detection of parallax), the problem of systematic residuals was not qualitatively ameliorated. This created something of a crisis. 
We realized that further improvements would be possible if we allowed for binary orbital motion. Lacking apparent alternatives, 
we went ahead and introduced projected binary orbital motion. This led to radical changes in our now (1 l+l)-parameterization. 
Allowing two dimensions of binary motion meant that both the orientation and the size of the binary separation could change. 
The latter induced changes in both the size and shape of the caustic, and so made it essentially impossible to define the time 
parameters At cc in such a way that was at the same time mathematically consistent and calculable in a reasonable amount of 
time. We therefore went back to something very like the original geometric parameterization but with four additional parameters 
(dt , q, a, uo , £e , , P* , fE , d, uf). Here the direction of 7Te contains the information on the direction of the binary axis relative to 
the line of ecliptic latitude at time t , when the source is closest to the binary center of mass in the Sun's frame of the reference. 
The binary orbital motion is incorporated via d and w. The quantities, a, u , and t^ are all given in the frame of the reference of 
the Sun as well. 

However, while this parameterization has the advantage of mathematical simplicity, it would have introduced severe insta- 
bilities into the fitting procedure. At bottom, the problem is the same as the one that led to the substitution (to, p*) — ► (t cc ,i, Ati) 
described above, but substantially more damaging. This is because the microlens parallax is a relatively poorly constrained quan- 
tity. Yet, within the framework of this parameterization, a change of the trial value of 7Te of only 1% (by itself) would shift t CC) i 
by ~ 0.3 day, and so induce A% 2 ~ 1000. Thus, to avoid such huge x 2 jumps, even more careful choreography would have been 
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required. 

Instead, we made the following changes to the parameter scheme. First, we made the reference time, t c , be the time of 
the closest approach of the source to the cusp, rather than the closest approach to the center of mass. Second, we adopted, 
for the frame of reference, the frame of the Earth at t c rather than the frame of the Sun. Our final choice of parameters is 
(dt c , q, a! ,u c ,t' E ,t c , p*, 7Te, d, uj), where u c is the impact parameter relative to the cusp and a 1 and t' E are evaluated at t c and in 
the frame of reference of the Earth at the time. The change of reference frame is responsible for the form of the parallax deviation 
given in equation (8), but from a practical point of view it is very helpful so that the non-parallax parameters do not change very 
much when the parallax is changed: in particular they are similar to the solution without parallax (i.e. 7Te = 0). The particular 
choice of the reference time t c is useful because most of the "action" of the event happens close to this time, either during the 
cusp approach itself, or during the second caustic crossing a week previously. Hence, both t c and u c are relatively well fixed by 
the data, while the angle a is also relatively well fixed since it is strongly constrained by the cusp approach and second crossing. 
As a consequence, we are able to find relatively robust minima for each (d, q) grid point in about a single day of computer time, 
which is quite adequate to reach a global minimum. 

Unfortunately, the (1 l+l)-parameterization failed to qualitatively lessen the problem of systematic residuals. We then rec- 
ognized that more LD freedom was required, and so added a square-root LD parameter in addition to the linear one. This reduced 
the systematic residuals to < 1%. As a result, we fit limb darkening with a two-parameter form so that the final fit we adopted is 
an (1 l+2)-parameterization fit. 



D. Hybrid Statistical Errors 

As we described in Appendix C, our x 2 minimization procedure is effective at fixed (d, q) (with the nine other geometrical 
parameters allowed to vary), but does not work when all 1 1 geometrical parameters are allowed to vary simultaneously. We 
therefore find the global minimum by evaluating x 2 over a (d, q) grid. How can the errors, and more generally the covariances, 
be determined under these circumstances? 

In what follows, the parameters will be collectively represented by a vector a, and the indices i, j will be allowed to vary 
over all p — 11 parameters. We assign a\ — d and a 2 = q, and designate that the indices to, n will be restricted to these two 
parameters. What we seek to evaluate is Cjj = cov(ai, aj) = (diCij) — (ai)(a,j). 

First we note that it is straightforward to determine c mn : simply evaluate x 2 at a series of points on the (ai, (12) grid, and fit 
these to x 2 = Xmin + El, n =i b mn (a m - a™ m )(a„ - < lm ). Then (c mn ) = (5 TO „) _1 is the covariance matrix restricted to the 
first two parameters, i.e., c mn — c mn . 

Next, at fixed (a\, a 2 ) = (a\, a5,), we evaluate the restricted covariance matrix of the remaining p — 2 = 9 parameters by 
varying 0, 1, or 2 parameters at a time and fitting the resulting x 2 hypersurface to x 2 — X m in + Sij=3 bij{ a i ~ a f m )( a i ~a™ m )- 
Then, using the result derived in the appendix of Gould & An (2002), one can show that (cy) = (bij)^ 1 is related to the full 
covariance matrix by 

Cij = Cjj ^ ' b mn c m iC n j 7^ 1, 2) , (Dl) 



rn.n 



and the parameters, di, at the constrained minimum are 



E« in - 4> m „c m (i^l.2). (D2) 

m.n 



Differentiating equation (D2) with respect to yields 



dat 



Eddi 

dan 

The partial derivatives can be determined simply by finding the change in as one steps along one axis of the (d, q) grid. Since 
c mn is already known from the first step, above, the c mi are also known. Finally, the remaining covariances can be found by 
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substituting equation (D4) into equation (Dl), 



^ mn n f a t ' 

oaL a a: 
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Table 1 . PLANET /-band Photometry of EROS BLG-2000-5 



Telescope 


Error Cut 
(mag) 


Seeing Cut 
(arcsec) 


Number of Points 


Error Scaling 


Median Seeing 
(arcsec) 


SAAO 


< 0.03 a 


< 2.1 


428 


1.99 


1.41 


YALO 


< 0.03 a 


< 2.3 


424 


1.61 


1.58 


Canopus b 






333 


2.96 


2.82 


Perth b 




< 3.1 


161 


3.63 


2.44 



a The formal value reported by DoPHOT. 

b The difference imaging analysis result has been used. 



Table 2. Relations between Parameterizations 



Physical Parameters Fit Parameters 



tc 


tc 


Q 


q 


p* 




Te 


+k,±) 1/2 




dtc^-E ("E, || &w - 7I"E,J_e n ) 




7r E 1 [( < ^ 7r E,|| + wdt c 7r E ,j.)e w + {-dir E . ± + wrf tc 7r Ei j_)e„] 
"cusp + u c 7TE 1 [(7r E) j.cosa / - 7r E) || sina')^ + (7r E) || cos a' + 7r Ei _L sin a')e n ] 


«tc 




t^T 1 7r E : 1 [(7r E j| cos a' + 7r Ei j_ sina')e w — (7t Ej _l cos a' — 7r E _|| sina')e n ] 



Note. — For simplicity, the reference times, t c , for both systems are chosen to be the same: the time of the closest approach 
to the cusp, i.e. (u tc — w C us P ) ■ u tc — 0. The additional transformation of the reference time requires the use of equations (6) 
and (8). Figure 2 illustrates the geometry used for the derivation of the transformation. The unit vector e„ points toward the NEP 
while e w is perpendicular to it and points to the west (the direction of decreasing ecliptic longitude). 
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Table 3. PLANET Model Parameters for EROS BLG-2000-5 







iinoprtaintv 3 ' 




1.928 


0.004 


q 


0.7485 


0.0066 


a' 


74? 18 


0?41 


u c 


-5.12 x 1CT 3 


3. x 10" 5 


t' E 


99.8 days 


1.5 day 


to 


1736.944 b 


0.005 day 


p* 


4.80 x 10~ 3 


4. x 10- 5 




-0.165 


0.042 




0.222 


0.031 


d 


0.203 yr- 1 


0.016 yr- 1 


LO 


0.006 rad yr" 1 


0.076 rad yr" 1 


ME.mj 


3.83 yr" 1 


0.49 yr" 1 


^E,n 


-2.82 yr" 1 


0.44 yr" 1 


ME 


4.76 yr- 1 


0.13 yr- 1 


a ec c 


-36? 3 


7? 8 


tte 


0.277 


0.008 


d 


-126?5 


3? 6 



a l-(T errorbar. The uncertainties of fit parameters are determined by fitting \ 2 distribution to a quadratic hypersurface. For 
more details, see Appendix D 

b Heliocentric Julian Date - 2450000. 

The angle of fi E with respect to ecliptic west. 

d The angle of d tc with respect to ecliptic west. 
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Table 4. Limb-Darkening Coefficients for EROS BLG-2000-5 



rv 


0.452 ±0.075 


A/ 


0.011 ±0.137 


Ti cos $ + A/ sin $ a,b 


0.207 ±0.156 


A/ cos $ — T/ sin $ c ' b 


0.402 ±0.003 


C/ 


0.552 ±0.090 


di 


0.011 ±0.139 


ci cos * + di sin * ad 


0.290 ±0.166 


d/ cos * - c/ sin * c - d 


0.470 ±0.003 



a rotational transformation that maximizes the variance 

b $ = -61?32 

c rotational transformation that minimizes the variance 

d * = -57?14 

Note. — The error bars account only for the uncertainty in the photometric parameters restricted to a fixed lens model, 
determined by the linear flux fit. 



Table 5. Kinematic Characteristics of the Lens and the Source 





location 




Ox 










(kms- 1 ) 


(kms" 1 ) 


(kms- 1 ) 


(kms- 1 ) 


Vh,p 


disk 


-18 


38 


25 


20 


v S,p 


disk 


-11.5 


31 


21 


17 




bulge 




100 


100 


100 



a asymmetric drift velocity 



Note. — The x-direction is toward the Galactic center from the LSR, the y-direction is the direction of the Galactic rotation, 
and the z-direction is toward the north Galactic pole. The lens is assumed to be an M dwarf while the source is a K3 giant. The 
quoted values for the disk components are derived from Binney & Merrifield (1998). 
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Fig. 1. — PLANET /-band light curve of EROS BLG-2000-5 (the 2000 season only). Only the data points used for the analysis 
("cleaned high-quality" subset; see § 2) are plotted. Data shown are from SAAO (red circles), Canopus (blue triangles), YALO 
(green squares), and Perth (magenta inverted triangles). All data points have been de-blended using the fit result - also accounting 
for the seeing correction- and transformed to the standard / magnitude; / = I B — 2.5 log[(F(t) — Fb) / F s ]. The calibrated source 
magnitude (I s = 16.602) is shown as a dotted line. The three bumps in the light curve and the corresponding positions relative 
to the microlens geometry are also indicated. 
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Fig. 2. — Geometry used for deriving the transformation shown in Table 2. The direction of dt c is chosen to be the x-axis while 
7Te lies parallel to the direction of decreasing ecliptic longitude; e w . The reference time, t c is the time of the closet approach to 
the cusp, M CUS p- 
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Fig. 3. — Prediction of deviations of light curves for SAAO and YALO from the geocentric light curve for a chosen model. The 
solid curve is the magnitude difference between the SAAO light curve and the geocentric one, and the dotted curve is the same 
for YALO. Nominal night portions (between 6 pm and 6 am local time) of the light curve are highlighted by overlaid dots. 
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Fig. 4. — Close-up of model light curves for the end of the second caustic crossing. The solid curve is modeled for SAAO 
observations, the dotted curve is for YALO, and the dashed curve is the geocentric light curve. The timing of the end of the 
second crossing for SAAO is earlier than for YALO by 1 1 minutes. For comparison purposes, all the light curve are calculated 
assuming no blend. 
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Fig. 5. — Surface brightness profile of the source star. The thick solid curve is the prediction indicated by the best fit model (see 
Table 4). In addition, the variation of profiles with the parameters allowed to deviate by 2-er along the direction of the principal 
conjugate is indicated by a shaded region. For comparison, also shown are theoretical profiles taken from Claret (2000). The 
stellar atmospheric model parameters for them are log g = 1.0, [Fe/H]= —0.3, and T c g — 3500 K (dotted curve), 4000 K (short 
dashed curve), 4500 K (long dashed curve). Note that the effective temperature of the source is reported to be 4500 ± 250 K by 
Albrow et al. (2001b) and 3800 ± 200 K by Minniti et al. (2002). 
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Fig. 6. — "Magnitude" residuals from PLANET model of EROS BLG-2000-5. Symbols are the same as in Fig. 1 The top panel 
shows the light curve corresponding to the time of observations, the middle panel shows the averaged residuals from 15-sequential 
observations, and the bottom panel shows the scatters of individual residual points. 
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Fig. 7. — Same as Fig. 6 but focuses mainly on the "anomalous" part of the light curve. The middle panel now shows the daily 
averages of residuals. 
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-1 1 

9/9 E 

Fig. 8. — Geometry of the event projected on the sky. Left is Galactic east, up is Galactic north. The origin is the center of mass 
of the binary lens. The trajectory of the source relative to the lens is shown as a thick solid curve while the short-dashed line 
shows the relative proper motion of the source seen from the Sun. The lengths of these trajectories correspond to the movement 
over six months between HJD = 2451670 and HJD = 2451850. The circle drawn with long dashes indicates the Einstein ring, 
and the curves within the circle are the caustics at two different times. The solid curve is at t — t c while the dotted curve is at 
the time of the first crossing. The corresponding locations of the two lens components are indicated by filled (t = t c ) and open 
(the first crossing) dots. The lower dots represent the more massive component of the binary. The ecliptic coordinate basis is also 
overlayed with the elliptical trajectory of 7Te<; over the year. 
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Fig. 9. — Close-up of Fig. 8 around the cusp approach. The source at the closest approach (t = t c ) is shown as a circle. The solid 
curve is the caustic at t = t c while the dotted curve is the caustic at the time of the first crossing. The positions of the source 
center at the time of each of the observations are also shown by symbols (same as in Figs. 1, 6, 7) that indicate the observatory. 
For the close-up panel, only those points that were excluded from the fit because of numerical problems in the magnification 
calculation (see §§ 2 and 4) are shown. Note that the residuals for all points (included these excluded ones) are shown in Figs. 6 
and 7. 




Fig. 10.— Calibrated CMD for the field around EROS BLG-2000-5, observed from YALO. The arrow shows the reddening 
vector. The positions of the source (S), blend (B), and the center of the red giant clump (C) are denoted by capital letters. 
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Fig. 11. — Kinematic likelihood for v as a function of Ds- The three curves are for different distributions of the source velocity: 
near-disk-like; v$ — v Iot + vs, p (short-dashed line), bulge-like; vg — us.p (solid line), and far-disk-like; v$ = —v Iot — v$. p 
(long-dashed line). The top panel shows the likelihood derived using only the two-dimensional projected velocity information 
while in the bottom panel, the likelihood also includes the radial velocity information derived from the high resolution spectra of 
Castro et al. (2001). 
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Fig. 12. — Distributions of the ratio of transverse potential energy, |Tj_| = [q/ (1 + q) 2 ]GM 2 /r± to transverse kinetic energy, 
K± = [q/(l + q) 2 ]Mv 2 _/2, for binaries seen at random times and random orientations, for three different eccentricities, e = 
0, 0.5, 0.9. Also shown is the 1 a allowed range for EROS BLG-2000-5. Non-circular orbits are favored. 



